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Abstract 

NASA Glenn Research Center, ONERA and the University of Illinois are conducting a 
major research program whose goal is to improve our understanding of the aerodynamic 
scaling of ice accretions on airfoils. The program when it is completed will result in 
validated scaled simulation methods that produce the essential aerodynamic features of the 
full-scale iced-airfoil. This research will provide some of the first, high-fidelity, full-scale, 
iced-airfoil aerodynamic data. An initial study classified ice accretions based on their 
aerodynamics into four types: roughness, streamwise ice, horn ice, and spanwise-ridge ice. 
Subscale testing using a NACA 23012 airfoil was performed in the NASA IRT and 
University of Illinois wind tunnel to better understand the aerodynamics of these ice types 
and to test various levels of ice simulation fidelity. These studies are briefly reviewed here 
and have been presented in more detail in other papers. Based on these results, full-scale 
testing at the ONERA FI tunnel using cast ice shapes obtained from molds taken in the IRT 
will provide full-scale iced airfoil data from full-scale ice accretions. Using these data as a 
baseline, the final step is to validate the simulation methods in scale in the Illinois wind 
tunnel. Computational ice accretion methods including LEWICE and ONICE have been 
used to guide the experiments and are briefly described and results shown. When full-scale 
and simulation aerodynamic results are available, these data will be used to further develop 
computational tools. Thus the purpose of the paper is to present an overview of the program 
and key results to date. 


I. Introduction 

In many applications, simulation of an ice accretion geometry on a wing or other surface is required for 
aerodynamic evaluation. The measurement of aerodynamic performance from an aircraft or wing with actual ice 
accretions is very difficult and expensive and seldom a practical solution. Aircraft performance with ice accreted in 
flight are available, 1 " 3 but limited, due to the cost and difficulty in acquiring these data. In flight it is difficult to 
document the ice accretion conditions and geometry accurately. For accretions formed in an icing wind tunnel the 
inability to provide adequate instrumentation, nonuniformity of the cloud over the model, poor aerodynamic flow 
quality, and cost also make aerodynamic measurements difficult and rare. The most common way to acquire iced- 
airfoil and wing data is to use simulated accretions in a dry-air wind tunnel. 
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The best current technology for accurate ice accretion simulation is the mold and casting method developed at 
NASA Glenn. 4 In this method molds are made from an ice accretion generated in an icing wind tunnel, such as the 
NASA Icing Research Tunnel (IRT). From these molds ice accretion castings are made that maintain the major 
features of the ice including the detailed surface roughness and the spanwise and chordwise variations. These 
castings are then typically attached to wings and airfoils, and instrumented to obtain high-fidelity aerodynamic 
data. 5 ' 7 This is an expensive process and not practical for many situations, but does generate benchmark data for 
iced-airfoil and wing research. 

While much of our understanding of ice accretion aerodynamics is anchored by cast ice simulation data, there 
have been very few studies employing full-scale ice accretion cast shapes aerodynamically tested at full-scale 
Reynolds number and Mach number. The second FAA/NASA Tailplane Icing Program provides full-scale, cast-ice 
performance data at full-scale Reynolds and Mach number, but the ice shapes were applied to a fully three- 
dimensional tailplane model. 5 The NASA Modern Airfoils Program 9 currently provides the best airfoil performance 
data set at flight Reynolds and Mach numbers, but due to size limitations in the NASA Low-Turbulence Pressure 
Tunnel (LTPT) used for the aerodynamic testing, the ice accretions were obtained on a 3-foot (0.91-m) chord airfoil. 
Ice accretion scaling methods can provide scaling of the gross ice shape for moderate scales, 10 but not of the 
roughness and ice shape details known to be aerodynamically important in some cases. 

The simulation of ice accretions on airfoils and wings is needed or desirable for many applications. For flight 
tests or wind-tunnel tests of iced airfoils for certification or research, cast ice shapes are seldom available and 
simulations are required. Often wind tunnel testing requires scale models for which no castings are available and 
small-scale simulations are needed. Computational methods need to model ice on airfoils and wings and it is either 
impossible or impractical to model a full, rough 3-D accretion. Simpler geometries that accurately represent the key 
aerodynamic features are needed. Currently we have only a limited understanding of how to accurately simulate the 
aerodynamic effect of ice accretion on lifting surfaces. Accurate simulations need to be based on an understanding 
of the fundamental flow physics for the simulations to be robust under a wide range of geometries and flow 
conditions. Here we describe a computational and experimental program using scale and full-scale aerodynamic 
data to develop the physical understanding and scaling methods for ice-shape classifications based on known flow 
physics. 

Thus the research program described in this paper addresses the iced-airfoil aerodynamic simulation and scaling 
issue. The program is a joint effort between NASA Glenn, the French research organization ONERA, and the 
University of Illinois. The goal is to provide for the first time, high-fidelity, full-scale, iced-airfoil aerodynamic data 
and validated scaled simulation methods that produce the essential full-scale aerodynamic features. Preliminary 
results of this program have been presented in previous papers. 1 114 The purpose of the current paper is to present an 
overview of the program and key results to date. 

To optimize our use of the expensive full-scale aerodynamic facility, ice accretion and aerodynamic testing 
were done in scale initially. The NACA 23012 airfoil was chosen for this research program as representative of a 
family of airfoils commonly used in industry and has shown in many cases to be aerodynamically sensitive to ice 
accretion. Thus aerodynamic penalties are amplified and a more severe testing of the simulation methods results. 


II. Initial Research Results and Plans 

A research program was designed to provide a validated ice accretion simulation method for iced-airfoil 
aerodynamics. The plan includes both full-scale and subscale ice accretion and aerodynamic wind-tunnel testing in 
several phases. By using subscale testing whenever possible to minimize the expense of full-scale testing, and as a 
way to validate the scaled ice accretion simulations, an efficient and cost effective program was designed. 

Phase I: Ice-Shape Classification 

The first phase of the research was to develop the ice-shape classifications based on the aerodynamics unique to 
each type. There are many types of ice accretion documented in the literature categorized based on the accretion 
process. These include rime ice, glaze ice, mixed ice, beak ice, runback ice, intercycle ice, etc. While these 
classifications are appropriate for thinking about the accretion of ice, they may not be as useful when the objective is 
aerodynamic simulation. In an initial stage of this research Bragg, Broeren and Blumenthal 12 examined our 
understanding of iced-airfoil aerodynamics and classified icing into four types based on the flowfield physics. 
Much of this discussion is taken from Bragg et al. 12 and more details on these classifications and iced-airfoil 
aerodynamics can be found there. Figure 1 provides a qualitative picture of the four classifications and their 
geometry and aerodynamic effect. 
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The first classification to be discussed, ice roughness, occurs during the initial stages of the accretion process 
before a significant ice shape, such as a horn, is accreted. The other three ice types are also “rough,” but here the 
initial surface roughness does not significantly alter the airfoil contour and thus the inviscid flowfield. Ice 
roughness may be associated with glaze, rime or mixed conditions and also from ice-protection system operation. A 
key aerodynamic feature of ice roughness is that it is usually much larger than the local boundary-layer thickness, at 
typical flight Reynolds number, even at the very early stages of development. These roughness elements then each 
act as bluff bodies with their own 3-D separated flowfield. In fact, the distinguishing characteristic of ice roughness 
and the other types of ice can be described in terms of the extent of boundary-layer separation. Ice shapes of the 
streamwise, horn or spanwise ridge categories are characterized by separated flow regions that are primarily 2-D 
with characteristic lengths that are large compared to the height of the ice accretion. In contrast, the separated flow 
regions generated by the ice roughness are fundamentally 
3-D, are very local to the roughness elements, and are of 
similar scale as the roughness itself. 

Roughness is characterized by its height, density (or 
concentration) and surface location. The effect on airfoil 
performance is dependent upon all of these parameters. 

Roughness shape can also be significant, but of the 
irregular shapes seen in ice roughness it is not thought to be 
as important, and is certainly less well understood than the 
other three parameters. Roughness affects airfoil 
performance by first directly increasing the skin friction. 

Roughness can also cause early boundary-layer transition 
and promote thickening of the boundary layer leading to 
early trailing-edge separation. These effects then manifest 
themselves through modified skin friction and pressure 
distributions into performance degradation — increased drag 
and decreased maximum lift. 

Longer icing exposures can cause roughness to build 
into horn ice, usually associated with glaze, or clear, icing 
conditions. The horn-ice shape can be characterized by its 
height (k). angle it makes with the chord line (9) and its location indicated by sic. the non-dimensional surface 
length. The most pervasive feature of the horn-ice flowfield is the large separation bubble that forms downstream of 
the horn. The stagnation point is usually located on the ice shape and the boundary layer cannot negotiate the large 
adverse pressure gradient encountered at the tip of the horn. Thus the separation location remains relatively fixed at 
the horn tip over a large angle of attack range. The separated shear layer undergoes transition to turbulent flow and 
then usually reattaches to the airfoil surface downstream. The presence of this separation bubble causes a large 
redistribution of pressure that results in pitching moment changes and decreased lift. The presence of the bubble 
greatly increases the airfoil drag. This flowfield is also known to be unsteady and three-dimensional. These aspects 
make computational simulation extremely challenging. Therefore, understanding the behavior of the separation 
bubble is key to understanding horn-ice aerodynamics. 

For this shape it is somewhat paradoxical that the separation bubble features that makes the flow complicated 
for analysis and computational simulations at the same time simplifies the effects on lift, drag and pitching moment. 
For example, several studies 15 " 17 have shown that very simple geometries, such as a leading-edge spoiler can 
reproduce the performance characteristics of a horn ice shape. In fact, this is addressed later in this paper and by 
Busch et al. 18 in much more detail. The point is that a simple geometry representing the height, angle and location 
of the ice hom, essentially generates an equivalent separation bubble on the airfoil, and hence, very similar 
performance results. Also, the addition of surface roughness is shown in some studies 6 ' 15 ' 18 ' 19 to have only minor 
effects on the integrated performance. 

Streamwise ice has the smallest effect of the non-roughness ice accretions as depicted in Fig. 1 and has received 
the least attention in the literature. Therefore, an understanding of the aerodynamics of these accretions is not as 
developed as for the horn and spanwise-ridge ice. Streamwise ice is often formed as a result of rime icing 
conditions that occur at cold temperatures when the incoming droplets freeze on the surface upon impingement. As 
a result of this phenomena, the initial ice accretion forms following the contour of the airfoil surface. At large 
accretion times, or when the icing conditions are otherwise appropriate, streamwise-ice shapes can occur that are not 
as conformal to the original airfoil surface and may grow a horn-like feature into the flow. In some cases flow 
separation may result as with horn ice. However, these separation bubbles tend to be much smaller and therefore 



Fig. 1 Qualitative description of aerodynamic 
effects for various ieed-airfoil flowfields. 12 
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have less of an effect on overall flowfield and aerodynamic performance. This means that other flowfield features, 
such as trailing-edge separation play at least an equal role in the aerodynamics. For streamwise ice, the specific ice 
geometry and surface roughness can be important factors in the aerodynamics. 

For the streamwise ice geometries that are conformal to the airfoil leading edge, the stagnation point at 
moderate lift coefficients occurs on the ice shape and the boundary layer remains attached as it flows around the 
leading edge of the ice and downstream on the upper surface. Since the streamwise ice/airfoil intersection is not 
smooth, an adverse pressure gradient may exist in this area and flow separation may occur in the junction region. 
The flow separation location is not fixed to a specific point on the ice shape as it is to the tip for a horn shape, but 
can move upstream or downstream depending upon the angle of attack and incoming boundary-layer state which 
depends on the surface roughness, Reynolds number, etc. For a streamwise ice shape having a horn-like feature 
oriented into the flow direction, a larger separated flow region may exist especially at higher angle of attack. These 
separation bubbles are typically much smaller than for horn ice and thus do not have as large of an effect on the stall 
mechanisms. 

Spanwise-ridge ice accretions are perhaps most often associated with SLD-type icing conditions. Usually these 
accretions form downstream of leading-edge ice-protection systems and, despite their association with SLD, can 
occur for all drop size ranges. Runback icing can form ridge accretions and occurs when there is a heated leading- 
edge ice -protection system not operating at 100% evaporation. Water flows back on the surface from the heated 
section to freeze downstream on the cooler, unheated surface. Due to the formation mechanisms just described, 
ridges often exhibit extensive spanwise variation in their geometry. These properties, and the associated flowfield, 
make the spanwise ridge-type accretion very different from the horn shapes discussed previously. As indicated in 
Fig. 1, the spanwise-ridge ice can have more severe aerodynamic effects and the shapes themselves are typically 
more 3-D than horn shapes. Spanwise ridges are generally located farther downstream than horns or streamwise 
shapes. This distance allows the boundary layer to develop and perhaps transition or to become transitional due to 
small ice roughness upstream of the spanwise ridge. Thus the spanwise ridge has the characteristic of a flow 
obstacle. 

The unique characteristics of the ridge, large spanwise variation and downstream location, make the resulting 
aerodynamics and performance different from the horn-ice case. The only similarity is the large separation bubble 
downstream of the spanwise ridge that can be complex, 3-D and highly unsteady. Spanwise ridges also have a 
separation bubble upstream of the ridge because the ridge is located well downstream of the stagnation point. This 
upstream separation can be a further challenge for computational simulation. The fact that the ridge is located well 
downstream of the stagnation point also amplifies the effect of airfoil geometry, characterized by a clean airfoil 
pressure distribution, much more so than in the horn case. 

The aerodynamics of spanwise ridges are dependent upon the ridge geometry, size, location and airfoil section. 
In the case of geometry, Lee and Bragg 20 tested four different shapes of simulated spanwise ridges: forward facing 
quarter round, backward facing quarter round, half round and forward facing ramp all having the same maximum 
height on the NACA 23012m airfoil. The result was a large variation in iced Ci max from 0.25 for the airfoil with the 
forward facing quarter round and ramp to up to 0.55 for the airfoil with the half round simulation. Perhaps a much 
more important result is the effect of location and airfoil geometry. Lee and Bragg 20 showed the variation in 
maximum lift penalty with the spanwise-ridge chordwise location and related the most critical location to the 
pressure distribution on the clean airfoil. 

This brief summary of the salient characteristics of each ice-shape type presents a convenient classification of 
iced-airfoil aerodynamics. Real ice shapes may not fit neatly into only one classification and of course, it is not 
realistic to think of these classifications as rigid, or unchanging. Undoubtedly more research will yield greater 
insight into the aerodynamics and may lead to further development of these concepts. Perhaps an additional 
classification will be added in the future. Future considerations notwithstanding, the present analysis formed the 
foundation for the research described in this paper. 

Phase II: Subseale Model Ice Accretion Testing 

The objective of this phase was to obtain high-fidelity ice shapes having characteristics of the four types 
described in Phase I. To obtain these ice shapes, a subscale model was subjected to simulated, in-flight icing 
conditions in an icing wind tunnel. The icing conditions were based upon the in-flight icing environment a 
commuter aircraft might encounter. Molds were made of the ice accreted on the model under these conditions. 
Castings were then made from these molds and used as the high-fidelity ice shapes for the tests in Phase III. Note 
that the purpose here was not to produce scaled ice accretion geometry, but to generate ice molds that were 
representative of the four ice classifications to use for aerodynamic simulation development in Phase III. 
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The subscale ice accretions were acquired in the NASA Glenn Research Center’s Icing Research Tunnel (IRT). 
The NACA 23012 airfoil model used in this phase of the testing had a 1. 5-foot (0.46-m) chord and spanned the 
height of the IRT test section (6-feet, 1.83-m). The model was equipped with pressure taps that facilitated model 
alignment and a removable leading edge section that facilitated the molding of selected ice accretions. For the 
spanwise-ridge ice accretion, an electric foil heater was attached to the removable leading edge. Icing conditions 
and heater settings were adjusted to create a situation where water would mnback over the heater and then freeze as 
a spanwise ridge downstream of the heater on both the upper and lower surfaces of the model. A photograph of the 
model in the IRT is shown in Fig. 2. 

Selecting the icing conditions for this phase of the 
program began by considering a full-scale commuter 
aircraft operating in various natural icing environments. 

A representative commuter aircraft with a wing section 
similar to the NACA 23012 airfoil was reviewed and 
analyzed in terms of aircraft weight, airspeeds, lift 
requirements, etc. The analysis provided information on 
the desired airspeeds and model angles of attack to use in 
the testing. Icing cloud conditions were initially selected 
from the FAA’s Appendix C icing envelopes. Icing 
scaling methods were then applied to help attain 
reasonable similarity in impingement limits, collection 
efficiency, freezing fraction, accumulation parameter, etc. 

Further adjustments to the icing conditions were made so 
that the testing conditions fell within the IRT’s operating 
envelope. 

Analysis of the testing conditions was conducted 
using the LEWICE ice accretion software program and is 
discussed in more detail later in Section VII. LEWICE cases were run for the various icing conditions for both the 
full-scale and subscale models. This analysis helped to ensure the resulting ice accretions were both a) 
representative of actual inflight icing conditions and b) exhibited the geometric features that would produce the four 
classifications of ice shapes determined in Phase I of this study. 


Table 1. Subscale Ice Accretions Selected for Aerodynamic Testing 


Ice 

Classification 

Run 

No. 

Airspeed 

knots 

AoA 

deg. 

MVD 

(.im 

LWC 

g/m 3 

Total 

Temp. 

F/C 

Static 

Temp. 

F/C 

Spray 

Time 

min. 

Roughness 

ED0762 

200 

2.0 

15.4 

0.75 

28.0/-2.2 

18.5/-7.5 

0.5 

Horn 

ED0735 

200 

2.0 

15.4 

0.75 

28.0/-2.2 

18.5/-7.5 

5.0 

Streamwise 

ED0730 

175 

5.3 

15 

0.30 

0.0/-17.8 

-7.4/-21.8 

5.0 

Spanwise Ridge 

ED0760 

175 

0.9 

15 

0.64 

24.0/-4.4 

16.8/-8.5 

5.0 



Fig. 2 Subseale NACA 23012 ice accretion model 
installation in NASA IRT. 


The ice accretions produced in the IRT were recorded using photographs, ice tracings, and ice depth 
measurements. When a particular ice accretion was selected for Phase III testing, the accretion was repeated and a 
mold was made for that run. Molds were made of a total of sixteen ice accretions. Six of these were horn 
accretions, four were streamwise accretions, and there were three each of both the roughness and spanwise ridge 
accretions. Castings were subsequently made from each of these molds. The best representatives of the four 
classifications of ice shapes determined in Phase I then formed the reference ice shapes for aerodynamic testing in 
Phase III. Table 1 provides the test conditions for these four ice accretions. 
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Phase III: Subscale Model Aerodynamic Testing 

Phase III of the program focused on exploring simulation 
methods for selected ice accretions in each of the four 
classifications. Since each of the four classifications has different 
fundamental aerodynamics, if techniques can be developed to 
simulate these four shapes, most ice accretions found in practice 
can be simulated. A variety of simulation methods can be found in 
the literature and in practice and were categorized for the subscale 
testing. The simulation categories or methods considered were: 3- 
D castings, 2-D smooth, simple geometry and simple geometry 
with spanwise variation. To account for the roughness associated 
with ice accretion, distributed roughness may be added to each of 
the non-casting simulation methods. The highest fidelity 
simulation method is considered to be the 3-D casting, since this is 
manufactured directly from a mold of the ice accretion. It is 
considered the benchmark for aerodynamic data that the other 
simulation methods are evaluated against. The purpose of this 
subscale testing was to explore many variations of the simulation methods in order to evaluate their ability to 
reproduce the proper iced-airfoil aerodynamics. In this way, the accuracy of the various simulation methods was 
systematically quantified. Quantifying this accuracy at small scale and understanding the aerodynamic differences 
of the simulation methods is key to successfully executing Phase VI of this program. 

The ice accretions 
selected in each of the four 
categories are shown in 
Figs. 3-6. Figure 3 shows 
the ice accretion typical of 
initial roughness (cf Phase 
II). There is a very smooth 
zone on the leading edge in 
the stagnation-point region, 
followed by distributed 
roughness downstream. 

Tracings for this type of ice 
shape were not obtained due 
to the inability of the ice 
tracing method to 
adequately capture such 

geometric features. The horn shape shown in Fig. 4 is a classic glaze ice shape with an upper and lower horn 
structure typical for this icing condition. The streamwise shape shown in Fig. 5 is typical of a rime ice condition, 
having a smooth zone in the stagnation region with rime feathers and roughness downstream. Finally, the spanwise- 
ridge shape is shown in Fig. 6. There is both an upper and lower-surface ridge that formed due to freezing of water 
running downstream off of the heated surface. Casting simulations were manufactured for aerodynamic testing from 
molds of each of these accretions. 

The various lower-fidelity simulation methods are described in detail by Busch et al. IS and are summarized here 
using the horn-ice shape as an example. The 2-D smooth simulation has been used in previous work. 6 ' 9 A digitized 
pencil tracing of the ice shape was processed using Smagglce." 1 This NASA software assists in removing point 
anomalies that often occur in digitized pencil tracings. It can also be used for smoothing, however, for this example 
no smoothing was performed with Smagglce. This tracing was then “extruded” across the model. The resulting 
shape is shown in Fig. 7. The simple-geometry simulation method entails replacing the major ice features with basic 
geometric shapes. An example of this for the horn example is shown in Fig. 7. The upper and lower surface horn 
structures have been replaced with rectangular shapes that mimic the height, angle and surface location. This 
simulation method is also based on previous work by Papadakis, 1517 who used “spoiler-ice” simulations and Kim 
and Bragg. 22 Other variations of the simple geometry were also considered including a simulation where the horns 
had a periodic spanwise variation in height based on measurements of the casting simulations. Both the 2-D smooth 
and the simple-geometry simulations were also tested with various sizes of roughness applied. 




Fig. 3 Photograph of roughness shape 
ED0762. 
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The subscale model 
aerodynamic testing was 
carried out at the University of 
Illinois using the low-speed, 
low-turbulence wind tunnel. 
This facility has a 3 -ft by 4-ft 
test section capable of speeds 
up to 160 mph. An 18-inch 
chord NACA 23012 model 
was designed and built 
specifically for this program. 
For this model, the Reynolds 
number was 1.8x1 0 6 , 

corresponding to a Mach 
number of 0.18. The lift and 
pitching moment were 

acquired from a floor mounted 
force balance. Surface static 
pressures were measured and 
also used to determine the lift 
and pitching moment. The 
drag was measured via wake 
survey using momentum- 
deficit methods. 

Figure 8 shows a picture 
of the NACA 23012 model 


Fig. 6 Tracing and photograph of spanwise ridge shape ED0760. 




installed in the University of 
Illinois wind tunnel. The 
model was designed with 
interchangeable leading edges 
that facilitated the various ice 
simulations. There was a 
baseline leading edge having 
the NACA 23012 profile that 
was used to document the 
uniced, or clean, airfoil 
performance. There were two 
ice leading edges having 
truncated nose geometry that 
allowed for the attachment of 
the ice simulations. The 
simulations bolted on to this 
leading edge and thus had a 
rigid, repeatable mounting 
system. The model had a main 
chordwise row of pressure 
taps, a secondary chordwise 

row and a set of spanwise taps on the upper surface. The ice-shape simulations also had pressure taps installed. In 
preparation for this work, a detailed study was made regarding pressure tap installation in ice simulations and those 
results were applied to this work. 13 14 More details of the simulations, experimental methods and model can be 
found in Busch et al. 18,23 



Fig. 7 Comparison of 2-D smooth and simple-geometry simulation for horn 
shape ED0735. 
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A summary of the aerodynamic results for the four 
ice shapes listed in Table 1 is presented in Fig. 9. The 
ice simulations in this case were all castings of the ice 
accretion, considered the highest fidelity. The data 
show the striking similarities in the integrated 
performance for the airfoil with roughness shape 
ED0762 and the streamwise shape ED0730 (cf. Figs. 3 
and 5), despite the differences in size, shape and icing 
conditions. This was likely due to the location of both 
ice shapes being on the airfoil leading edge. At this 
location, the size of any ice roughness (or small 
accretion) is large compared to the local boundary-layer 
thickness. Kim and Bragg 22 observed similar results 
for geometric ice-shape simulations located on the 
leading edge of an NLF-0414 airfoil. The horn shape 
had the largest effect on C/ and C m , particularly evident 
in the C„, was the deviation from the clean values 
resulting from the large redistribution of surface 
pressure by the separation bubble(s) aft of the horn(s). Fig. 8 Subscale NACA 23012 aerodynamic model 

It is interesting that while the horn shape had the largest installation in University of Illinois wind tunnel, 

effect on C/ and C m , the drag penalty was less than for 

the spanwise-ridge shape over the interval -4 < a < 4 deg. Outside of this interval, the spanwise-ridge shape had a 
lower overall drag penalty than for the airfoil with the horn shape. The C m characteristics for the spanwise-ridge 
shape illustrate the effect of the larger size ridge located on the lower surface that resulted in much greater deviation 
in C„, at negative angle of attack. The stall characteristics for the airfoil with this shape are also more ambiguous 
than for the other shapes and are reminiscent of the thin-airfoil stall type. 24 Judging by the break in the C m curve, the 
stall angle appears to be near 8 deg., which was accompanied by an attendant slope change in the C/ curve, while C/ 
continued to increase with angle of attack up to 14.4 deg. 

The aerodynamic effect of various fidelity simulations is shown for the streamwise-ice shape in Fig. 10. In 
terms of lift variation with angle of attack, similar results were observed for the castings and simple-geometry 
simulations. The airfoil with the 2-D smooth simulation had a slightly higher C tmax . The simple-geometry 
simulation yielded pitching moment behavior that was closer to the uniced airfoil than for the airfoil with the casting 
simulation, while the 2-D smooth simulation more closely matched the casting. This same trend was also reflected 
in the drag data over the interval -5 < a < 9 deg. These data indicate that the ice roughness associated with the 
casting simulation, but not present on the other two simulations plays an important role in the aerodynamics. This is 
consistent with the notion that large-scale boundary-layer separation is not a significant flowfield feature for this ice- 
shape classification. Instead, characteristics like trailing-edge separation are likely more significant. In this case, 
the amount of roughness present can affect the boundary layer downstream of the ice simulation leading to different 
trailing-edge separation behavior. This effect is shown in Fig. 11 where there is a deviation of upper-surface 
pressure for the airfoil with the casting simulation from the clean airfoil aft of x/c = 0.60. This deviation in pressure 
likely indicates boundary-layer separation. 

The effect of the addition of roughness to the 2-D smooth simulation on the integrated performance is illustrated 
in Fig. 12. Glass microbeads with k/c = 0.0009 (35-40 grit size) applied to the 2-D simulation resulted in very 
similar iced-airfoil performance to the airfoil with the casting simulation over a large angle of attack range. 
Increasing this roughness size to k/c = 0.0033 and 0.0026 (14 and 20 grit carborundum grains, respectively), resulted 
in more severe performance penalties. What is not clear, and is a part of this ongoing research, is what scale sets the 
size of the roughness. That is, looking at the roughness on the casting, it appears to be of a physical size and shape 
more akin to the k/c = 0.0033 and 0.0026 carborundum grains. The k/c = 0.0009 microbeads seem much smoother 
than the ice roughness on a qualitative level, yet produce the best match in terms of integrated performance. 
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Fig. 9 Aerodynamic performance comparison for the NACA 23012 model with the casting simulations. 
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Fig. 10 Aerodynamic performance comparison of simulation fidelity for streamwise ice shape ED0730 on the 
NACA 23012 model. 
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It is possible that there are also roughness concentration effects that may be significant, as well as more complicated 
coupling between the roughness and gross ice shape effects. 

Performance comparisons for initial ice 
accretion roughness simulation are 
summarized in Fig. 13. In this case there is ' 

no 2-1) smooth simulation because the ice , [ j j J a Casting 

accretion was simply classified as I i | □ 2 -d smooth 

roughness. So the various simulations were _ 40 ; \ 1 | I j I I i 

k/c = 0.0009 glass microbeads, and k/c - C \ \ I 1 I 

0.0033 and 0.0026 carborundum grains. In f 3 0 i/uv< J j j j 

the simulations, the smooth zone near the - 

stagnation point (cf. Fig. 3) was preserved _ 20 :.! \\\ | J [ j j j 

and the roughness was applied to the upper : ; I j 

and lower surface. The lift and pitching , ,, I | j j 

moment data indicate some sensitivity to : j, 

the roughness size, with the k/c = 0.0026 00 ; 1 I i — j — | - -Ss 

case having the best overall match. The : J 

real challenge for roughness simulation, Tl. I i J j j I 

however, lies in the drag. As the data o'o' ' ' 'o.V ' ' ’oV ’ ' 'o! 3 ' ' ' o!7 ' ' o'.s ' ' ' 0.6 ' ' 'oV ' ' oV ' ' ' 0.9 Vo 

indicate, it was difficult to match the drag x/c 

increase resulting from the casting 

simulation on the airfoil. The drag data Fig. u Pressure distribution comparison for the casting and 2- 
indicate that the k/c = 0.0009 roughness D smooth streamwise ice shape ED0730 simulations on the 

simulation offers the closest match to the NACA 23012 model at a = 11.4 deg. 

casting for lower angle of attack while the 

k/c = 0.0026 roughness simulation offers the closest match at higher angle of attack. These data suggest that further 
experiments be conducted where different grit sizes are used on the upper and lower surface. In this case, using k/c 
= 0.0026 roughness on the upper surface and k/c = 0.0009 on the lower surface seems like the first logical 
combination. Furthermore, estimates of the actual roughness size on the casting indicate that this offers the best 
geometric size match as well. 


Fig. 11 Pressure distribution comparison for the casting and 2- 
D smooth streamwise ice shape ED0730 simulations on the 
NACA 23012 model at a = 11.4 deg. 



a (deg) a (deg) 


Fig. 12 Aerodymanie performance comparison of roughness applied to the 2-D smooth streamwise ice shape 
ED0730 simulations on the NACA 23012 model. 
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a (deg) a (deg) 

Fig. 13 Aerodynamic performance comparison for various roughness shape ED0762 simulations on the 


NACA 23012 model. 



a (deg) a (deg) 

Fig. 14 Aerodynamic performance comparison of simulation fidelity for spanwise-ridge ice shape ED0760 on 
the NACA 23012 model. 

The simulation results for the spanwise-ridge category of ice shape are summarized in Fig. 14. The lift and 
pitching moment data indicate that the 2-D smooth simulation applied to the airfoil resulted in a much better match 
than the simple-geometry simulation to the results from the casting simulation applied to the airfoil. This supports 
the notion of Bragg et al. 11 ' 12 that the specific cross-sectional geometry of a ridge simulation is important as regards 
its affect on the airfoil aerodynamics. As with many of these cases, drag coefficient simulation was much more 
challenging. The 2-D smooth simulation tested on the airfoil resulted in lower drag on the interval -2 < a <9 deg. 
than for the airfoil with the casting simulation. Outside of this interval the opposite was true. Attempts to improve 
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the simulation accuracy with the addition of roughness are illustrated in Fig. 15. It is interesting to note that it was 
the largest roughness size applied to the 2-D smooth simulation that resulted in the best matching of C,j, while the 
smaller roughness sizes resulted in higher drag penalties. At the same time, the addition of this roughness reduced 
the agreement in lift and pitching moment. Further, similar counter-intuitive results were observed in the variation 
of maximum lift with roughness size applied to the 2-D smooth simulation. Namely, smaller roughness sizes 
resulted in larger C/ max penalties. Clearly, these data illustrate the simulation challenges associated with the 
spanwise-ridge category of ice shape and how common expectations can be inaccurate. 


• — 

— Clean 

A 

— Casting 

B 

— 2-D + k/c = 0.0009 roughness 

— e — 

— 2-D + k/c = 0.0026 roughness 

© 

— 2-D + k/c = 0.0033 roughness 




a (deg) a (deg) 

Fig. 15 Aerodymanic performance comparison of roughness applied to the 2-D smooth spanwise-ridge ice 
shape ED0760 simulation on the NACA 23012 model. 


The aerodynamic data for the various levels of ice-shape simulation fidelity gathered in Phase III aided in 
confirming and further refining the ice-shape classifications described in Phase I. The data illustrate the range in 
performance data for typical simulation methods as applied to a subscale model. Establishing the relationship 
between simulation fidelity and the resulting aerodynamic performance is a key goal of this program. 
Accomplishing this goal for the subscale model led to increased understanding of the effects and were directly 
applied to the full-scale model research phases. This is described in Phases IV and V. 


Phase IV: Full-scale Model Ice Accretion Testing 

The objective of this phase of the program was to obtain molds of ice accretions of each of the four ice 
classifications on a full-scale airfoil. In a scaling methods development effort such as this, it is important to have a 
set of benchmark data and firm knowledge of the full-scale phenomena. This information provides validation data 
for both the subscale testing and scaling methods developed. In order to gain full-scale aerodynamic effects 
information for the four classifications of ice, accurate representations of the ice needed to be obtained. Therefore, a 
full-scale ice accretion model was built and tested in the IRT. This effort provided the full-scale, reference ice 
shapes for the program. 

The model built for this phase was again a two-dimensional, NACA 23012 that spanned from floor to ceiling of 
the test section of the IRT. It was machined from solid aluminum and had a removable leading edge section for 
making molds of the ice accretions. However, it had a chord length of 6.0-feet (1.83-m); four times that of the 
subscale model. A photograph of the model mounted vertically in the IRT is shown in Fig. 16. It was instrumented 
with pressure taps for alignment purposes and with thermocouples to monitor model temperature. Similar to the 
subscale model, an electric foil heater was attached to the removable leading edge. Icing conditions and heater 
setting were adjusted to create a situation where water would runback over the heater and then freeze as a spanwise 
ridge downstream of the heater on both the upper and lower surfaces of the model. 
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Test conditions for the full-scale icing tests were selected as indicated 
in the Phase II discussion; i.e., icing conditions that a commuter aircraft 
might encounter in flight as defined by the FAA’s Appendix C icing 
conditions and that would yield ice accretions falling within the four ice 
classifications described in Phase I. Some adjustments in the test 
conditions were made such that they fell within the IRT’s operating 
envelope. The full-scale icing conditions are given in Table 2. 

The full-scale ice accretions produced in the IRT were, as in the 
subscale testing, recorded using photographs, ice tracings, and ice depth 
measurements. When a particular ice accretion was selected for the Phase 
V testing, the accretion run was repeated and the removable leading edge 
and resulting ice accretion were detached from the model. The leading 
edge and ice were installed in a mold box and the mold was poured. Molds 
were made of a total of ten ice accretions. There were three each of the 
horn and streamwise accretions and two each of the roughness and 
spanwise ridge accretions. Castings were again subsequently made from 
each of these molds. The castings were then prioritized based upon which 
were the best representatives of the four ice classifications and then upon 
which were of most interest for the purposes of forming a full scale, high- 
fidelity, ice contaminated aerodynamic validation database. These ice 
shapes and test plans are described in the following section. 


Table 2. Icing Conditions for Full-Scale Icing Tests 


Ice 

Classification 

Run 

No. 

Airspeed 

knots 

AoA 

deg. 

MVD 

(ini 

LWC 

g/m 3 

Total 

Temp. 

F/C 

Static 

Temp. 

F/C 

Spray 

Time 

min. 

Horn 

EG 1 164 

175 

5.0 

20 

0.85 

28.0/-2.2 

20.8/-6.2 

11.3 

Streamwise 

EG1162 

150 

2.0 

30 

0.55 

-8.0/-22.2 

-13.5/-25.3 

10.0 

Roughness 

EG1126 

200 

2.0 

20 

0.50 

28.0/-2.2 

18.6/-7.4 

2.0 

Spanwise Ridge 

EG1159 

150 

1.5 

20 

0.81 

20.0/-6.7 

14.8/-9.6 

15.0 

Streamwise 

EG1125 

200 

2.0 

15 

0.30 

4.0/-15.6 

-5.3/-20.7 

20.0 

Roughness 

EG1134 

200 

2.0 

40 

0.55 

4.0/-15.6 

-5.3/-20.7 

2.0 



Fig. 16 Full-scale NACA 23012 
ice accretion model installation in 
NASA IRT. 


Phase V: Full-scale Model Aerodynamic Testing 

In Phase V of the program, aerodynamic testing will provide the icing community performance data for an 
airfoil with fill-scale, high-fidelity, ice accretion simulations. As was the case for the subscale testing, ice shapes 
were selected according to the four classifications developed in Phase I. In contrast to the subscale testing, only the 
high-fidelity ice casting simulations will be tested on the full-scale model. The reasoning here is that this approach 
increases the number of different ice shapes that may be tested with the given amount of tunnel time. A total of six 
ice accretions were selected for aerodynamic testing and are shown in Figs. 17-22. The first four of these shapes 
(Figs. 17-20) are from each of the four ice-shape classifications and were accreted using the icing conditions given 
in Table II. The two additional shapes (Figs. 21 and 22) were selected in part because they can be challenging from 
the standpoint of simulation on a subscale model. 

The aerodynamic testing will be carried out in the ONERA FI wind tunnel shown in Fig. 23. This is a large 
pressure tunnel having a test-section 3.5 m high by 4.5 m wide. The 6-foot (1.83 m) chord NACA 23012 airfoil will 
be oriented vertically and supported by a force balance located beneath the floor. Like the subscale model, it has 
two interchangeable leading edges. The baseline leading edge has the clean NACA 23012 profile and the ice 
leading edge is truncated such that the ice-shape simulations can be bolted in place on the nose of the model. Lift 
and pitching moment data will be acquired from the force balance as well as from the integration of the surface 
static pressures. Drag will be obtained via wake survey. The pressure tunnel will be fully utilized to explore any 
Reynolds number variations (from 4.0xl0 6 to 12xl0 6 and from 8.0xl0 6 to 16.0xl0 6 ) at fixed Mach numbers (0.10 
and 0.20, respectively), as well as Mach number variation (0. 10 to 0.28) at fixed Reynolds number (12.0xl0 6 ). 
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The full-scale model 
testing will be comprised of 
two campaigns. The first 
campaign is currently 
scheduled for spring of 2007 
and the second for summer 
2007. The objective in the 
first campaign is to measure 
the aerodynamic performance 
(lift, drag, pitching moment, 
surface pressures, etc.) for 
each of the six ice shapes over 
a large angle of attack range 
including stall. Flow 

visualization will also be 
performed. Fluorescent mini- 
tufts will be attached to the 
model upper surface over a 
segment of the span away from 
the pressure measurement and 
wake survey stations. The 
mini-tuft response will be 
recorded on video during the 
angle of attack sweeps. 

Surface oil flow visualization 
will be performed for selected 
cases as time allows. The 
objective of the second 
campaign is to acquire detailed 
flowfield data for 

computational methods 

development and validation. 

In this case, the time- 
dependent velocity field in the vicinity of an ice shape simulation will be acquired using an optical method such as 
particle-image velocimetry (PIV). Because of the time consuming nature of these measurements, only one ice shape 
will be selected for investigation. 

Phase VI: Simulation Validation Testing 

The full-scale model aerodynamic testing conducted in Phase V 
will provide a substantial validation data set for the final phases of 
this program. In Phase VI, the ice shapes used in the full-scale test 
will be scaled and simulated for testing on the subscale model using 
the Illinois low-speed facility. This will effectively “close the loop” 
by quantifying the accuracy of 2-D and simple geometry simulations 
for subscale iced-airfoil aerodynamic testing at low-Reynolds 
number. 

The current testing plan is to fabricate and test lower fidelity 
simulations on the subscale NACA 23012 model of the ice shapes 
shown in Figs. 17 to 22. The aerodynamic data, including flow 
visualization will be used to quantify the uncertainty with which full- 
scale iced-airfoil effects can be performed on a subscale model. This 
will provide full-scale to subscale comparison of iced-airfoil 
aerodynamics under tightly controlled conditions. 



Fig. 19 Picture of roughness shape 


EG1126. 






- 0.05 0.00 0.05 0.10 0.15 

x/c 


Fig. 17 Tracing and photograph of horn shape EG1164. 
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Fig. 18 Tracing and photograph of streamwise shape EG1162. 
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Phase VII: CFD Research 

In addition to the 
experimental effort currently 
underway, computational studies 
are also being performed. The 
objectives of this effort are to use 
current CFD capabilities to assist 
the experimental effort and to 
use data from the project to aid 
in validation of 3D CFD methods 
as applied to ice-shape 
geometries. Additionally, it is 
desired to use the experimental 
data from the project to evaluate 
the requirements for geometric 
simulation by ice-accretion 
codes. The latter objective is 
based upon identifying the 
geometric characteristics that 
influence the aerodynamics and 
using that information to 
determine what elements of an 
ice-accretion simulation must be 
accurately recreated by such a 
computational method. 

In order to accomplish these 
objectives, a series of 
computational tasks has been 
identified. These consist of two- 
dimensional simulation of the ice 
shapes generated in the NASA 
Icing Research Tunnel (IRT), 
and two- and three-dimensional 
CFD calculations of the 
aerodynamics of the measured 

ice shapes. The ice-accretion simulation task is designed to compare 
the results of the computation with the ice shapes produced in the 
NASA IRT. Additionally, the ice-accretion simulation tool was used 
to aid in selection of the icing conditions for the test program. Two- 
dimensional CFD calculations of the ice shapes generated from the 
previous task as well as for the respective ice-shape tracings from the 
experiment will be performed. These calculations will be compared 
to understand the differences between ice-shape tracings and 
computer generated ice shapes with respect to the resulting 
aerodynamic forces. These results will also be used to evaluate the 
aerodynamic differences between sub-scale and full-scale ice shape 
geometries. Three-dimensional CFD calculations of selected ice- 
shape/airfoil geometries that have been tested in the experiments will 
be performed. The results of these calculations will be compared to 
the data provided from the experiment. Analysis of the 
computational results will be performed to determine what elements of the calculation have reproduced the 
experiment and what deficiencies, if any, have been identified. 



Fig. 22 Picture of roughness shape 
EG1134. 
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Fig. 20 Tracing and photograph of spanwise-ridge shape EG1159. 
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Fig. 21 Tracing and photograph of streamwise shape EG1125. 
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Ice accretion simulation results 

Both the NASA LEWICE code 25 and the ONERA 
ONICE code 26 have been used to calculate ice 
accretions in support of this program The LEWICE 
code was used to simulate the ice shapes produced for 
the program during testing in the IRT. This code 
produces two-dimensional ice shape profdes using 
input airfoil geometry, flight conditions, and icing 
cloud conditions. The code calculates ice shapes by 
performing a flowfield analysis, calculations of water 
droplet trajectories, a mass and energy balance at the 
icing surface, and an ice growth calculation based upon 
the mass of ice accumulated and its distribution on the 
airfoil surface. The ONERA code uses a finite element 
method with a C grid to solve the velocity potential 
equation to determine the flowfield and with 
Lagrangian droplet tracking to determine trajectories 
and impingement locations. The Messinger model is used for ice-accretion thermodynamics and correlations are 
used for convective and evaporative heat transfer. The Makkonen correlations enable a representation of roughness 
effects. The ice shape is first estimated for a given icing time in one step. Then the flowfield, the trajectories and the 
heat transfer coefficient are calculated for this ‘estimated shape’. Assuming that the values of the local collection 
efficiency and heat transfer coefficients vary linearly from their values on the clean airfoil to their values on the 
profile covered by the estimated shape, the thermodynamic balance is made and the final ice shape is calculated. 

The initial use of the LEWICE code was to aid in selection of the icing conditions that were used in the 
experimental testing in the IRT. As mentioned previously, the code was used to help determine the icing conditions 
that would produce ice shapes that exhibited the geometric features of two of the four classifications of ice shapes 
determined in Phase I of this study. LEWICE was not used to produce surface roughness or the ice ridge shapes 
which were produced experimentally in the IRT with the use of an electrical heater blanket on the leading edge to 
produce runback. 

In conjunction with the actual IRT testing, LEWICE calculations were performed to compare the code results 
with the corresponding measured ice shapes. Typical results are shown in Figs. 24 and 25 which show comparisons 
for the 18- inch chord model and 72-inch chord model, respectively. Since Table 2 presents accretions where molds 



Fig. 23 ONERA FI wind-tunnel test section. 



Fig. 24 LEWICE and ONERA simulations of IRT test using 72-inch chord airfoil model. 
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a) horn case (ED0735) 



Fig. 25 LEWICE and ONERA simulations of IRT test using 18 inch chord airfoil model. 


were taken and therefore no tracings, calculations were performed for similar runs where tracings of the shapes were 
available. Thus in Fig. 24 case a) the horn case (run EG1163 which is similar to 1164 except MVD=20, LWC=.85 
and spray time=l 1.2) and b) the Streamwise case (run EG1 161 which is a repeat of 1162). In Fig. 25 the conditions 
are those of Table 1 with a) the Florn case (run ED0735) and b) the Streamwise case (run ED0730). The Span 
notation in the legend refers to the spanwise location where the experimental ice shape measurement was taken with 
36-inch being the tunnel centerline. 

The results for streamwise shapes show excellent agreement and the results for horn shaped geometries show 
agreement typical for ONICE and LEWICE. In the latter case, this type of agreement has been shown 27 to be within 
the experimental repeatability of the ice shape generation process. Comparisons for roughness shapes were not 
performed because experimental tracings were not taken and the codes do not produce ice shapes with small scale 
roughness features. 

Aerodynamic analysis of ice shape profiles with CFD 

The data from this test program will ultimately be used, aside from the aerodynamic scaling analysis, for 
comparison with CFD results in order to assess the current computational capabilities as well as to aid in 
improvement of those methods. The CFD effort incorporates the use of both two-dimensional and three-dimensional 
computational methods. The two-dimensional computations will be undertaken first in order to provide initial 
evaluations of the IRT ice accretions. The performance of these two-dimensional models can then be evaluated with 
the data from the tests at Illinois and at the ONERA FI tunnel. Once a reasonable understanding of the two- 
dimensional results has been developed, an extension to three-dimensional computations can be undertaken building 
upon the two-dimensional analysis. 

Some initial calculations for the two-dimensional data have been performed using a commercially available grid 
generation tool, -8 GRIDGEN, and a publicly available CFD code, 29 ARC2D. Modifications to the ARC2D code 
make it attractive in the evaluation of flowfields about ice shape profiles. 30 The modifications were to the turbulence 
model to account for the well documented issues surrounding the Baldwin-Lomax model near regions of flow 
separation. In addition, the Elsa code 31 ' 32 was used by ONERA to perform iced-airfoil calculations with a k-omega 
turbulence model. This code provides the numerical simulation of the compressible viscous mono-species, for 
steady and unsteady flows, on three-dimensional or two-dimensional multi-block structured meshes. It is based on a 
finite volume approach. 

Some initial results are shown in Figs. 26 - 29 which show the grid and Mach number contours for clean and 
iced geometries. The selected geometries are those obtained by the LEWICE and ONICE codes. Note that run 
EG1 109 is a glaze condition similar to EG1 164 except LWC=0.5, T 0 =-6.7C, and spray time= 22.5. 

Future work in this project will include more two-dimensional aerodynamic calculations for the LEWICE and 
ONICE code generated shapes as well as the experimental ice shape geometries described previously. Additionally, 
more two- and three-dimensional calculations will be performed to identify the differences between two- and three- 
dimensional results. 
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Fig. 26 GRIDGEN Grid structures for (a) clean NACA 23012 airfoil geometry and (b) with a 
LEWICE generated ice shape. (EG1109) 




Fig. 27 


(a) (b) 

ARC2D Mach number contours for (a) clean NACA 23012 airfoil geometry and (b) with a 
LEWICE generated ice shape.(EG1109); a = 0°, Re = 3.29xl0 6 , M = 0.4. 
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Fig. 28 



x/c 

a) b) 

Grid structures (generated with ICEM) for (a) clean NACA 23012 airfoil geometry and (b) with an 

ONICE generated ice shape (Fig.24a). 



x/c x/c 

a) b) 

Fig. 29 Mach number contours (using the Elsa code) for (a) clean NACA 23012 airfoil geometry and (b) 
with an ONICE generated ice shape (Fig. 24 a); a = 0°, Re = 8.1xl0 6 , M = 0.3. 


III. Conclusions 

This paper provides an overview of the NASA/ONERA/lllinois research program to better understand ice 
accretion simulation for airfoil aerodynamics. Ice accretions have been classified into four groups by examining 
their aerodynamic effect: roughness, horn, streamwise, and spanwise-ridge ice. This program uses subscale and full- 
scale testing to provide dill-scale data and a development and validation of simulation methods of varying fidelity 
for these four classifications. Ice accretions were acquired in the NASA IRT on an 18-inch NACA 23012 airfoil. 
Guided by computational results, ice-accretion conditions were set to be representative of a commuter aircraft in 
Appendix C conditions. Molds were made for all four ice-shape types. Aerodynamic measurements in the Illinois 
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subsonic tunnel compared the performance and flowfield on the cast ice shapes to simulations of varying degree of 
fidelity. The height and location of the horn or ridge were dominate features on the horn and spanwise-ridge 
accretions. Here roughness played a secondary effect. Roughness was critical on the streamwise shape to reproduce 
the lift and the drag. The spanwise ridge was the most challenging to accurately simulate. Research is continuing in 
the Illinois tunnel to continue to improve our understanding of these effects. 

These small-scale simulation experiments are guiding the full-scale experiments. Ice accretions have been 
acquired on a full-scale NACA 23012 model in the IRT and aerodynamic measurements on cast ice shapes will 
occur later this year in the ONERA FI tunnel - a large pressure tunnel. Based on these results we will close the 
loop by using the knowledge gained to simulate these aerodynamic results in the Illinois tunnel using simulated 
small-scale shapes. This test will determine our ability to duplicate full-scale ice accretion in both geometric and 
aerodynamic scale (Reynolds number). In addition to providing the baseline for scaling development, the full-scale 
results will be used for CFD development and validation and contribute to a goal of identifying key ice shape 
characteristics and simulation methods for aerodynamic testing. 
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